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Abstract 

We study meson-baryon scattering with strangeness -1 in unitary chiral perturbation theory. Ten coupled 
channels are considered in our work, namely tt^A, tt^E", tt^S^ , tt"^ 

K^E~ . A large amount of experimental data are analyzed, including the recent precise measurement by the 
SIDDHARTA Collaboration of the energy shift and width of the Is state of kaonic hydrogen. This leads to a 
strong constraint on the free parameters in our theory and of the resulting meson-baryon scattering amplitudes. 
We also analyze the uncertainty that stems by using several different strategies to perform the fits to data. It 
is found that large uncertainties in the subthreshold extrapolation of the K~p scattering amplitude arise by 
employing either only one common weak pseudoscalar decay constant or distinguishing between /tt, /k and 
/,,. However, in both cases a good reproduction of experimental data is obtained. We also discuss the pole 
content of the resulting S-wa,ve amplitudes, particularly in connection with the two-pole structure of the A(1405) 
resonance. 
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1 Introduction 



Antikaon-nucleon reaction is quite an interesting subject in hadron physics. At low energies Chiral Perturbation 
Theory (ChPT) [U [21 13] constrains strongly the possible interactions but at the same time strong nonperturbative 
effects take place. In this respect we have the presence of the A(1405) resonance between the ttE and KN thresholds, 
which is also a manifestation of strongly coupled channel dynamics. There are also big SU{3) symmetry breaking 
effects, giving rise to a large split in thresholds because of the different masses of the coupled channels composed by 
the lightest octets of pseudoscalars and baryons. Nevertheless, in the chiral limit all these channels are degenerate 
which implies that despite the difference in the thresholds could be large all of them are relevant in a wide range 
of energies. By the same reason a realistic study of the resonances in meson-baryon scattering with strangeness —1 
in the energy region between 1 and 2 GeV requires to take all of them into account. 

In recent years, great progress by exploiting chiral effective field theory has been made on this research field 
[lEiaiTllHllSlIinillllllllllllllllSiniinillS- These studies corroborate the nature of the A(1405) [TO] as a 
dynamically generated resonance from KN and ttE interactions. However, in Ref. [5] it was already noticed the 
presence of indeed two nearby poles in the A(1405) region. This fact concerning the two-pole structure of the 
A(1405) has been studied in subsequent papers [3[8l|9l[l0l[Il[l3l[I3[l6l[IZl[20]. The two poles a,re commonly 
characterized as following: the first pole in the complex energy plane is located quite close to the KN threshold 
with a small imaginary part, around 10-30 MeV, and a strong coupling to KN. In turn, the second one is wider, 
with a relatively large imaginary part around 50-200 MeV, coupling more strongly to the ttS channel and its pole 
position shows more dependence on the specific theoretical model [101 [12] • 

On the other hand, the K~p scattering length can be determined by the measurement of the energy shift and 
width of the Is state in kaonic hydrogen. Due to the controversial results of the DEAR Collaboration [21] and the 
low precision of the KEK measurement [22] there was an uncertainty of around a factor 2 in this scattering length 
[7] [HI [H] • In this respect the recent precise measurement by the SIDDHARTA Collaboration [23] may finally fix the 
K~p scattering length, up to a precision of around 20%. Of course, an important point now is to reproduce this 
measurement simultaneously with all the other scattering data. Studies in this direction are already accomplished 
in Refs. [TS] [16] [TT] [18] . In the present work we give a step forward and take into account all the data included in 
Refs. [m [ini [17] [T5] and also consider additional data in our fits, such as the 7r~S"'" event distribution [2^, the 
cross section of K~p — >■ rjA [2S], measurements from the reaction of K^p — > 7r°7r°E° [20] and the it A phase shift at 
the S~ mass [271 [2H]) St^a- Moreover, to further constrain the fits we include the pion-nucleon isospin even S'-wave 
scattering length aQ_^, the nucleon a term a-^N and the masses of the N, A, S and S. We do not aim to give a 
precise description for a^^, (j-^n and the masses, since we only calculate them up to 0{p^) at tree level in SU{'i) 
baryon ChPT. Potentially large corrections from loops and higher order low energy constants could exist. So in the 
fit we allow relatively large errors for those quantities in order to account for our theoretical uncertainties. Even 
so, we still find that they can lead to important constraints on the free parameters. In addition, we discuss here 
two sources of ambiguity overlooked in previous studies. One point concerns with using either just one common 
weak decay constant for all the channels or distinguishing between /jr, Jk and The other controversial issue 
that also requires a closer look is the definition employed for the , as it has become customary in many cases to 
weight differently the data fitted corresponding to different observables. The first point is of relevance because it 
shows that the uncertainty affecting the subthreshold extrapolation of the K~p elastic S'-wave amplitude is much 
larger than estimated in the recent studies [151 HI] HZl HI] that also reproduced the SIDDHARTA measurement. 
Let us recall that this subthreshold extrapolation is of great interest for the active and controversial field on K 
few- nucleon systems and possible existence of bound states of K with heavy nuclei [211 1301 13IJ 132] • 

As an output of our study, we can also obtain important information on the rich resonant content or spectroscopy 
associated to strangeness —1 meson-baryon S'-wave scattering. In this respect, for J = we confirm the two-pole 
nature for the A(1405) resonance and find poles associated with the resonance A(1670). In / = 1 we find poles 
that are more dependent on the details of the fits. Poles in the region of the E(1620) are found. Other poles 
corresponding to the E(1750) and to more controversial 1 = 1 resonances around the KN threshold are obtained 
in one of the fits. 

After this introduction, we present the basic set up of our approach in Sec.[2J Next, we discuss in Sec.[3]the fits 
to the data included with the interacting kernel fixed at leading and next-to-leading order in the chiral expansion. 
We find a good reproduction of scattering data together with the energy shift and width of the Is state of kaonic 
hydrogen. Finally, we discuss in Sec. U] the pole content of our solutions. Section [5] includes a summary of the work 
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and our conclusions. 



2 Formalism 

Our starting point is the SU{3) meson-baryon chiral Lagrangian at 0{p) and 0{p^) |33) : 

£2 = bo{BB){x+) +bD{B{x+,B}) +bF{B[x+,B]+bi{B[u^,[u^\B]])) 

+ b2{B{Uf,,{u^,B}})+b3{B{u^,[u^,B]})+bi{BB){Uf,u^) + --- , (1) 

where we do not show those terms, represented by the eUipsis, that do not contribute to the 5'-wave meson-baryon 
scattering at tree leveL The octet of baryons N, A, S and S is cohected in the unitary 3x3 matrix B and the 
octet of hght pseudoscalar mesons tt, K, and ij is introduced through the basic chiral building blocks: the covariant 
derivative D^, x+ and u^. In our notation, — iu'^d^Uu^ U = = e*^*/-'' with $ the unitary 3x3 matrix for 
the octet of light pseudoscalars it,K, and 77. In this way the parameter / corresponds to the pseudoscalar weak 
decay constant / in the SU{3) chiral limit ^ij- Furthermore, Dfj^B = d^B + [F^, B] with F^ = [u^d^u + ud^u^)/2 
and x+ — '2.Bq{u^ MqU^ + uMqu) with Mq = {TO„,mrf,TOs} the diagonal mass matrix for the three light flavor 
quarks. Bq is defined in terms of the vacuum quark condensate in the SU{3) chiral limit as {0\q''q^\0) — —f'^B^S^^ 
[M] , The 0{p^) low energy constants bo, bo, bp, 6^=1^2,3,4 will be fitted to data. The values for D and F are taken 
from the determination of hyperon and neutron beta decays [3S]: D = 0.80 and F — 0.46. For more details on the 
notation and derivation of these chiral Lagrangians the reader is referred to [B^ . 

The perturbative amplitudes Vij for the meson-baryon scattering {(t>B)i — ^ {4>B)j up to 0{p^) have been given 
in Rcfs. [S HI m]. Here the labels i, j ^ 1,2,3, ... , 10 correspond to the ten different channels 7r°A (1), 7r°I]° (2), 
7r-E+ (3), TT+E- (4), K-^p (5), K"n (6), ryA (7), ryS^ (8), K°E° (9) and K+E- (10). Taking the initial center of 
mass (CM) three-momentum along the z-axis, the S'-wave amplitudes 7ij{W) can be projected out through 

%,{W)^ ^ J dnv,,{w,n,a,a), (2) 

where W is the energy in CM, a is the third component of the spin for the initial and final baryons and Q is the 
solid angle of the scattered final three-momentum. 

In SU{3) baryon ChPT the resummation of the unitarity chiral loops is crucial because of the large mass of 
the s quark and the corresponding larger masses of those pseudoscalar mesons and baryons with strangeness. 
Then, in many kinetic configurations the masses of the two particles in a given channel are much larger than the 
typical three-momentum and this enhances the two- particle reducible loop contributions [36]. In addition, we are 
also interested in the chiral dynamics involving the resonance region where the unitarity upper bound in partial 
waves could be saturated. Thus, it is not appropriate to treat unitarity in a perturbative way as in plain ChPT. 
The presence of the A(1405), which is near below the KN threshold, clearly signals that the perturbative chiral 
amplitudes can not be appropriate here and a nonperturbative method must be developed to probe the meson- 
baryon dynamics in the strangeness —1 sector. As a result, one must resum the right-hand cut that stems from 
unitarity and for that we employ Unitary ChPT (UChPT). This is based on an approximate algebraic solution 
to the N/D method [371 [351 132] , which was first applied to meson-meson interactions and then to meson-baryon 
scattering [Slllllli]. 

The unitarized meson-baryon scattering amplitude in this formalism can be cast in matrix notation as [B] 

T{W) ^ [1 + N{W) ■ g{W^)y^ ■ N{W) . (3) 

The function g{s), s — W^, collects the unitarity cuts contributed by the two-particle intermediate states and 
N{W) only contains the crossed-channel cuts. The unitarity loop g{s) can be calculated through a once subtracted 
dispersion relation (as well as in the dimensional regularization replacing the divergence by a constant) and for the 
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ith channel its explicit form reads 



16TT^g{s)i = aiip.) + log ^ - x+ log 
s + Mf - mf , 1 



1 x^ 
— X- log — 



(4) 



x± 



2s 



2\2 



± -^/-4s(Af2 _ ,0+) + (,s + Mf - m\) 



where ai (/i) is the subtraction constant that will be fitted to data, Mi and nii denote the baryon and meson masses 
in the ith channel. We point out that (/(s) is independent of the scale /x, introduced for dimensional reasons. 
Note that the combination ai(/Lt) + log mf//i^ could be reabsorbed in a scale independent subtraction constant. In 
our following discussions, we set /i — 770 MeV, that fixes the scale at which the subtraction constants ai(/^) are 
determined. We also constrain these subtraction constants by taking into account that when isospin symmetry 
is conserved the different attached to the states with different charges but made up from the same type of 
pseudoscalars and baryons should be the same. This was demonstrated for SU (3) symmetry in Ref. [TU] and we 
use its particularization for the isospin SU{2) subgroup. In this way, the three subtraction constants for 7r°S°, 
TT^E"*" and tt+S" are the same and equal to 02. The two subtraction constants for K^p and K'^n are denoted by 
05 and those for K^E^ and K^E~ correspond to ag. 

Concerning N{W) in Eq. ([3]), it does not contain any unitarity cut and in our present discussion it is simply 
equal to the perturbative partial wave amplitude T{W) in Eq. ([2]) calculated up to 0{p^). However, for a higher 
order calculation of the input chiral amplitude T{W), the interaction kernel N{W) is different from T{W) and the 
corresponding formula relating both can be found in Ref. [6]. 

With the previous set up, we are ready to calculate the cross section — >■ {<j)B)j, which in our normalization 
reads 

1 \Po\\rr 12 



(5) 



In the previous equation T(^,pB)i^(4>B)j corresponds to the unitarized S'-wave amplitude in Eq. Q, pi and pj 
denote the initial and final CM three-momentum of the baryons, respectively. We fit the cross sections of eight 



'P 



-Y.-), a{K-p 7r-E+), a{K-p 



different processes: a{K p ^ K p), a{K p — ^ K^n), (j{K 
a{K-p TT°A), a{K-p j]A) and a{K'p tt^tt^Y.'^). 

Two event distributions are also fitted: the 7r~E+ event distribution from K^p — > E+(1660)7r" 
and the tt^E'^ event distribution from K~p tt^tt^E" [55]. For the K~p — > tt'^tt'^E" reaction we consider the 
proton pole dominant exchange mechanism developed in Ref. [13], and evaluated as in Ref. [5], see its Fig. 30 
We normalize the tt'^TP event distribution by multiplying our theoretical results by a constant fixed in order to 
reproduce the highest peak value in the event distribution. 

To describe the 7r^E+ event distribution in the A(1405) energy region, we follow the scheme developed in 
Ref. [5] based on the assumption that the process takes places from an original scalar isoscalar source (with the 
same quantum numbers as the A(1405)) which is then corrected by final state interactions involving the three ttE 
and the two KN channels. For more details we refer to Ref. [Bj. Two extra free parameters r and r' are then 
introduced and the resulting formula is [B] 



dW 



32 



D- 



133 



D- 



134 



36, 



S+l 



(6) 



In this equation we indicate by D the matrix elements of the inverse of the matrix D=[l + T{W) ■ g(M^2)][2| 
Three ratios between different cross sections at the K^p threshold are also included in the fit [401 HTj : 



a{K^p — > 7^^E^ 
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Rr = 



(t{K p — > charged particles) 
(j{K-p all) 



Rn — 



a{K-p tt"A) 



a{K p all neutral states) 



(7) 



^In this exchange the proton is almost on-shell and this is the reason why the proton pole exchange mechanism is expected to be 
dominant 1131 . 

^Notice that the T-matrix Eq. JSj can be written as T(W) = D{W)-^N(W). 
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The formulas for ct^at, Qq^ and the masses of N, A, S and S at 0{p^) level calculated from the Lagrangians Ci 
and £2 in Eq. ([1]), read 





= -2 




2 


rriN 




rriA 


= mo 


ms 


= ma 




= Too 



(2&0 + 6d + &f) - {bi +b2 + 63 + 264) + 



mo - 2(60 + 2bF)M^ ~ 4(&o + bD- bF)Ml , 
/r2 



\,r2 



[bo - 2hF)Ml - 4(60 + bD + bF)Ml , (8) 

where we do not consider isospin-breaking effects for these quantities. The expressions for the masses of the baryons 
depend on the baryon mass in the SU{3) chiral limit denoted by too, which always enter as toq — 26o- Thus, from 
the study of masses one can not disentangle both parameters too and bo, although the latter can be fitted from 
scattering. In our study we fix toq = 0.9 ± 0.2 GeV from the higher order study [42], as in Ref. [8]. Its uncertainty 
is propagated to the errors in the final results obtained. Let us note that for the fits employing the interacting 
kernel calculated at 0{p) all the next-to-leading (NLO) counterterms 60, bo, bF and bi, i = 1,2, 3, 4, are set to zero 
in Eq. ^ and none of these quantities are then introduced in the fits at 0{p) described below. 

Finally, the new SIDDHARTA measurement on the energy shift AE and width T of the kaonic hydrogen Is state 
are related to the K^p scatting length aj^-p, after taking into account isospin-breaking corrections |43j . through 



AE- 



2o? fil aji-p [1 + 2a p,r (1 - In a) a^-p] 



(9) 



where a is the fine structure constant and fir — mpAIx- / {Mx- + mp) is the reduced mass of the K and p system. 
In our normalization the relation between the K^p scatting length and the unitarized amplitude Eq. ([3]) reads 

TK-p^K-piVs) 



Uk- 



(10) 



3 Fits and results 

Once we have introduced the theoretical formalism and the data included in the fits, we now proceed with the 
phenomenological discussions. In order to fix the free parameters we consider two different strategies to perform 
the fits, taking into account the fact that in literature there are two common ways to treat the light meson decay 
constants. For example, only one single decay constant is used for tt, K and 77 in Refs [H [SJ [SI El HI [131 HH] , while 
physical values are adopted in Refs. [HI [151 [HI [Hj. In Fit I, we use a common decay constant for all the states 
in the meson-baryon scattering amplitudes and then fit this single decay constant. In Fit II, we distinguish the 
decay constants for tt, K and 77, according to the pseudoscalar appearing in the state, so that one then employs 
Jk and frj, in order. The values = 92.2 MeV and Jk = HO.O MeV are taken from the PDG [45], while 

= 120 MeV (in fact it is J,,,) is estimated by using the result from Ref. [46]. Note that both strategies are 
compatible with the calculation of T in ChPT up to 0{p^), since the differences induced by using one or other 
decay constants are at least of 0{p^). For the masses of mesons and baryons, we always take their physical values 
in the scattering amplitudes in order to properly account for the threshold eflFects. 

To equally weight experimental data from different measurements, we divide the from each measurement 
by its number of data points and then sum over all the partial together. Then, the pcr degree of freedom 
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iXd.o.f) is defined as 



2 ^ L,fc"fc V- Xk 



( fh Gxp 



2 



X^=E^— ^F-^- (11) 

Here K is the number of different measurements considered in our fits, Uk stands for the number of data points in the 
fcth measurement, Up corresponds to the number of free parameters and y'f^l {yk%) represents the «th experimental 
(theoretical) point of the fcth set of data with standard deviation ak-i- This definition of is also used in 
Refs. [m im [ini im- This is basically done in order to enhance the weight of the kaon hydrogen data in the 
fits performed, because one assumes that this is a high quality data point. Nonetheless, we consider that using the 

/ 



previous definition of o / i'^ ^1- d instead of the standard one for a total set of X^a^i data points, 



Xd.o.f = Y] Xk : (12) 



k=l 

is somewhat arbitrary. In this way, we will also discuss the stability of the fits by switching from Eq. (jlip to 
Eq. dm). 

Concerning the experimental data, the references for the cross sections and event distributions are explained in 
Fig[TJ while the remaining were already introduced above. For the ratio 7 measured at the K~p threshold in Eq. (I?]), 
we notice that Ref. [40] reported the value 2.34 ± 0.08 and Ref. gj gave two values: 2.38 ± 0.04 and 2.35 ± 0.07. 
Other value 7 = 2.15 can be also found in Ref. [37]. Hence in our current work we decide to assign 7 a conservative 
error band at 5% level, i.e. 7 = 2.36±0.12, instead of 7 = 2.36±0.04 used in Refs. [HiadZKIl] and 7 = 2.38±0.04 
in Ref. [18 . A similar conservative error bar is also assigned for Rc, namely Rc = 0.664 ± 0.033. While since the 
figure Rn = 0.189 ± 0.015 derived in Ref. [3H] already contains an error bar larger than our conservative 5% level, 
we take this number with its error in our later fits|f| As mentioned previously, in order to account for our theoretical 
uncertainties, we assign large errors for the pion nucleon isospin even S'-wave scattering length a^^_, nucleon sigma 
term (Jt^n, and the masses of N, A, E and 3, whose central values and ascribed error bars can be found in the second 
column of Table [2] For ct^at one expects that it receives sizable higher-order corrections from the mesonic cloud 
which are expected to be positive and around 10 MeV [49^. On the other hand, the precise value for this observable 
is still under debate and one has the classic result a-^N = 45 ± 8 MeV [15] and the one favored from more modern 
experimental ttN data bases UTriv = 59 ± 7 MeV [SD]. This adds another 10 MeV of uncertainty, so that in our fits 
employing an 0{p^) calculation for this quantity we take cr^^v = 30 ± 20 MeV. Regarding the isoscalar scalar ttN 
scattering length its latest recent determination from pionic atoms is a,|_|_ = (7.6 ± 3.1) • 10"'^ 51 . Similarly 

as in Ref. [8] we expect a theoretical error in our 0{p^) calculation of around +1 • 10^^ n^^^i estimated from the 
0{p^) unitarity corrections [55]. As a result we take in our fits the value Qq^ = (0± 1) • 10^^ For the masses 

we include an error of 30% by considering a general expectation for the breaking of an SU{3) prediction. 

Now we are ready to present our fit results. The error bars of the parameters in the fits represent only the 
statistical uncertainty at the level of one standard deviation. For the definition of the in Eq. (fTT|) we take the 

distribution function for a system with i' degrees of freedom (being i/ the number of free parameters in the fit). 
For the other definition Eq. (jl2p . with a large numbers of degrees of freedom, we employ instead 

(13) 



as deduced in Ref. [53^ • Here, is the number of standard deviations, Xo is the minimum x^ calculated with the 
central values of the fit and results by taking a new configuration of free parameters. Those configurations with 
a x^ larger than the upper limit from Eq. (|13p with Ua- — 1 are rejected at the level of one standard deviation. The 
same set of parameter configurations will be also used in later discussions for determining the uncertainty in any 
observable considered. 




^We do not consider that we could have the phenomenology well under control within a precision of 5% or better. 
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3.1 Results using the defined in Eq. ( ITT]) 

We first discuss the results employing the definition for the corresponding to Eq. (|TT|) , in which every set of data 
corresponding to one observable is considered as an effective data point. The fitted free parameters are given in 
Table[T]for the different fits considered, namely, for the leading order fit, C'(p)-Fit, and for the NLO ones, Fit I and 
Fit II. From this table we observe that the hi coefficients have natural values which are quite similar between the 
two fits I and II (of course for the C'(p)-Fit they are fixed to zero). The largest differences happen for the values 
of the subtraction constants, specially for oi. Nonetheless, one has to stress that all the subtraction constants in 
our fits have natural size of 0(1). The cross sections and event distributions are shown by the (red) solid lines 
in Figs. [H [5] and [3] for the fits I, II and C'(p)-Fit, in order. The hatched area around each line is the estimated 
statistical uncertainty at the level of one sigma as explained above. In Table Owe give the reproduction for the rest 
of observables fitted, and show separately the energy shift and width of the Is state of kaonic hydrogen state in the 
left panel of Fig. H] by the empty symbols, as indicated in the figure. The rest of contents in Table [5] are shown by 
the empty squares and circles in Fig. [5] for fits I and II, respectively. In this figure every observable is appropriately 
rescaled so that all of them have similar size and the quality of its reproduction can be easily appreciated. 

As it is clear from Figs. [TJ [5]and Table[2]the reproduction of data for fits I and II is very good with &Xd o f lower 
than 1, see the first line in Table [TJ The reproduction of the observables is then consistent between each other. In 
this respect we find that our results are compatible with the value of S^^a from Ref. [55] and not with Ref. |27j . 
It is also remarkable that the 0(p)-Fit, having 6 free parameters less than Fit II and 7 less than Fit I, is able to 
achieve a good reproduction of data, with the important exception of the K^p — > rjA cross section that is not able 
to reproduce adequately. Indeed, this is the reason why the Xd o f ^he 0{p)-Fit is near 2. If the data of the cross 
section for K~p — )■ rjA is removed the resulting Xd o f 1.23, a much lower value. It follows then that the 0{p)-Fit 
can properly reproduce the lower-energy scattering data as well as the SIDDHARTA measurement on the energy 
shift and width of the Is state of kaonic hydrogen. As a result, one finds an indication towards convergence in the 
reproduction of all these data since it can already be well reproduced by the 0{p)-Fit (except the K~p — ^ ijA cross 
section) and the description is improved when including 0{p^) contributions to the interacting kernel (in which 
case the K~p — >■ rjA can also be accounted for). This indication of convergence in the reproduction of data when 
passing from the leading to NLO fits did not result when fitting the DEAR data [21] on the energy width and shift 
of the Is state of kaonic hydrogen |5]. We also find that the ratios in Eq. ([7|) are improved in Table [5] compared 
to the results obtained in Ref. f5^. The masses of the lightest baryon octet and (Jt^n are satisfactorily reproduced 
as well. Our results for AE and F in Eq. ^ are compatible with the SIDDHARTA measurements within errors, 
though our fits prefer somewhat lager central values for F, as shown in Table [3] and Fig. (T) It is clear from the 
right panel of Fig. [4] that the results from other recent studies [T5j [161 [ITj [18] also prefer larger values for F. 

We also calculate the K~p scattering length a^-p, Eq. (|10p . and those with / = and I — 1 extracted from 
data and the theoretical model. The results are shown in Table [3| Our values are compatible with the recent 
determinations from Refs. [1^1 [H]. In addition, we show in Fig. [5] the K~p — >■ K~p 5- wave amplitude around the 
threshold region and below it, the left panel is for the real part and the right one for the imaginary part. The points 
with error bars in Fig. [6| correspond to a^-p extracted directly from SIDDHARTA data by making use of Eq. ((9]), 
flif-p — (—0.65 ± 0.10) + i (0.81 ± 0.15) fm. We observe that this amplitude for fits I and II is quite the same at 
threshold and above it in the range shown in the figure, indeed the threshold parameters from both fits perfectly 
agree with each other in Table [3| However, quite different behaviors result in the energy region below threshold. 
The real parts of the K~p — )■ K~p scattering amplitudes from the two fits become incompatible below 1.42 GeV 
and the imaginary parts are also incompatible below 1.41 GeV. As a result, it is clear that a precise knowledge of 
the aj^-p scattering length, as that resulting from the SIDDHARTA data, can not pin down in a precise manner 
the extrapolation of the K~p 5-wave amplitude below threshold. Indeed, the difference from Fit I and II for the 
subthreshold extrapolation of the K~p S'-wave amplitude in Fig. [5] is much larger than the uncertainty estimated 
in Refs. [HI [13 [13 [TH] . However, an 0{p^) determination of the interacting kernel N{W) in UChPT, Eq. ©, 
could reduce this uncertainty because then one is sensitive to the change of the weak decay constants so that these 
changes have to be compensated by new terms at 0{p^). This is certainly an interesting calculation in order to 
improve the accuracy of our knowledge in this field of research. 
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Parameters 


Fit I 


Fit II 


0(p)-Fit 




X /d.o.j = 0.85 


xio.f =0-96 


xio.f = 1-87 


f (MeV) 


124.60ti'g8 


Fixed 


116.05+;-?? 

^^^ — 1.59 


bn (GeV^i) 


— 0.026 


-0.292+;"!!!? 
^ 0.007 





bn (GeV^^I 


-0.027+°-°?^ 

' —0.020 


o.ioi+!!-«J,« 

0.008 





bp (GeV~^) 


-,00+0.094 
^ — 0.082 


-0.200+°-?," 
0.011 





61 (GeV^^) 


714+0-»" 

' ^^-0.019 


n c:n9+0.005 





62 (GeV"^) 


1 001 +0.051 

1.00X__Q 


1 015+?-??t 

i-ui J_Q Q23 





b% (GeV"^) 




-0 306+°-°15 





64 (GeV"-'-) 


n Qon+0.037 
-U.»»y_Q 


„ onn+0.010 

-0.899+0012 





0.1 


-0.944 


. „„., -1-0.491 
4-761Io.397 


„ „77+1.509 
-<a.6{ (_i.288 


"2 


n 0^0+0.134 
-U.»ciO_(, 


-0.447_o,i63 


1 779+0.236 
^- ' '^-0.186 


CI5 


, n7Q+0.034 

l.U(0_Q QgQ 


1 fi«c:+0.042 
-i.b8t)_o 041 


-1.668t°;°S 


a? 




1 401+°-l**5 
i.1Ul_o 155 


r, 91 C- + 0.180 
Z.Z10_o i4g 


18 


1 QQQ+0.831 
i.aoo_j^ 4Q(j 


-0.1bS_o 050 


-0 170+°-2''^ 

U.i/U_o 239 


ag 


„ +0.035 
z.loi_Q (J22 


-2 406+°-°^** 

z.'lUD_o 027 


9 990+0.101 
z.zzci_o 068 


r (GeV-i) 


24.28t*;^« 


18.27i^;J« 


n 20+'*-^*' 


r' (GeV^i) 


10.85l^g;f2 


l'-00_14.06 


5 40+'^-^** 
■^•^^-18.94 



Table 1: Parameters from the three fits with the defined in Ea. (fTT|) . The -way in -which the error bars of 
parameters are calculated is explained in the text. 



Observable 


Input 


Fit I 


Fit II 


C>(p)-Fit 


AE (eV) 


283±36 


299I41 


2761^^ 




r (eV) 


541±92 


019+66 


60811° 


606^72 


7 


2.36±0.12 


2 36+°-" 
z.ou_o 22 


2 36+°-2'l 
z.oo_o 23 


2 36+°-2^ 

z.oo_o 28 


Rc 


0.664±0.033 


662+°-°°** 

U.DDZ_o.olO 


n CC1 +0.012 

0-661_o.oii 


647+0-006 
U.D4(_o 007 


Rn 


0.189±0.015 


1 09+0-025 

u.iyz_o.o2o 


n 1«S+0-028 
U.i8»_o.o29 


n iaQ+0.033 
U.i«S_o.o32 


5„A (degrees) 


3.2 ±5.3 




9 7+0.5 
' -0.4 


1 S+0-3 
--'--°-0.4 


a+_ {10-2m-i) 


Oil.O 


0.08t°-23 


46+°-22 
-U.4D_o.22 


66+°-°2 


a^N (MeV) 


30±20 




26.0t^°, 





Mn (GeV) 


0.94±0.28 


1.00t°-.f2 


Q2+°-20 

u.y^_o.2o 


mo 


Ma (GeV) 


1.12±0.34 


1 17+0.21 
-^-•-^'-0.21 


1 07+°-2° 

^-^'-0.20 


mo 


Ms (GeV) 


1.19±0.36 


1 14+0.20 
-^-•-^^-0.20 


1 19+0.20 
^-^^-0.20 


mo 


Ma (GeV) 


1.32±0.40 


1 oQ+0.22 
l.oo_o 24 


1 90+0.20 


mo 



Table 2: Results by using the fits in Table [TJ In the second column we give the values for the several observables 
as included in the fit. 



Observable 


Fit I 


Fit II 


0(p)-Fit 


Or-p (fm) 


-0 67+°-" + iO 92+°-°^ 
'-'-°'-0.12 + "-'■^^-0.08 


-o.6i«-.i? + ^ 0.891°;;;^ 


-0.61^°-.l|+^0.89^°-l 


0.1=0 (fm) 


_1 74+0.20 .,27+0.14 
-^--'^-0.17 * -^--^'-0.12 


-1.58t°-.^3«+^1.27lg-°^ 


-1.62l°-ji+il.32l°;i5 


oi=i (fm) 


30+012 , .„ .g+0.09 

'J-^^-o.io + ' '^-^"-o.og 


0.36t°-°^ + .0.52«-°^ 


40+°-°^ , -0 4^+0.10 



Table 3: K p scattering length, flj^-p, and / = 0, 1 KN scattering lengths, a/=o, ai=i, in order, for Fit I, Fit II and 
C'(p)-Fit given in Table [U 
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Parameters 


Fit I 


Fit II 


0(p)-Fit 




X /d.o.i = 1.87 


xio.f = 1-93 


xio.f =2-71 


f (MeV) 


^ 1 -^ — 0,94 


Fixed 


112.75+?-?^ 

1 ^ — 2.14 


bn (GeV^i) 


0.018 


-0.340+0-??:! 

^^^ — 0.009 





bn (GeV^^I 


— 0.012 


0.196+?-?" 





bp (GeV~^) 


-0.183+°-?,?? 

— 0.051 


-0.265+?-?;? 

0.014 





61 (GeV^^) 


856+°-?!;^ 


667+?-??! 

-0.007 





62 (GeV"^) 


1 621+0-?,?? 


1 193+?-?;? 

-^■-^^"^-0.017 





b% (GeV"^) 


-0 5oi+!;-?i? 


-0 504+?-?;? 

^•^^^—0.013 





64 (GeV"-'-) 


-, -, „„+0.026 

-1-177+0027 


1 non-t-0.012 
-1.082+0017 





0.1 


-0.800l|5;^g 


r)«i -1-0.391 
o.^41_Q 


c noQ+0.668 

-b.uad^o 5(jg 


"2 




n -1-0.103 

-U.00±_o 1Q2 


1 Kt;c; + 0.139 

-1.655_,o 135 


CI5 


-0 961+"-"^** 


1 5';8+o-o''i 

i.OOO_o 


-1.610lH?^ 


a? 


n 6^1+0-21'^ 

U.DOi_g j^yg 


, r,7c:-l-0.213 
l.Z(0_o 259 


-2.068t0;;» 


18 


-5.250_2 go7 


-0.168_o 


U.^lMO^o 221 


ag 


-2 160+°-"l2 


9 9QS+0-014 
z.zjo_Q 012 


2 165+0-0'''^ 

Z.1UO_0 QgQ 


r (GeV-i) 




19.78i^;«^ 


10.79l^5;fg 


r' (GeV^i) 




2i-76l^:i^ 


7 no+3.80 
' ■'^^-4.08 



Table 4: Parameters from the three fits with the defined in Ea. ([T2)) . The error bars are calculated as explained 
in the text. 



Observable 


Input 


Fit I 


Fit II 


C'(p)-Fit 


A-B (eV) 


283±36 


3341^1 


3241^1 


3081^^ 


F (eV) 


541±92 


649t«« 




664±fo 


7 


2.36±0.12 


r, r,i -1-0.40 

^•^-"--0.30 


9 17-1-0.37 

^•J-'~0.37 


2.00tO-.^3l 


Rc 


0.664±0.033 


652+0-007 


n 649+0 009 

U.D4y_0 oil 


640+0-00'' 
U.D4U_o 006 


Rn 


0.189±0.015 


p. 997-1-0.038 
u.zzi _o.o39 


n lQS+0-041 
U.19S_o.o29 


n 207+0-052 
^■^'J' -0.051 


SnA (degrees) 


3.2 ±5.3 


-i.ot;.? 


1 0-1-0.4 
^-•■^-0.4 


-1 7+0.3 


(10-2m-i) 


Oil.O 


O.lOtO-Ji 


1 oq-l-0.16 
--l-.oo_0.16 


-U. 'U_o.03 


CT.^]V (McV) 


30±20 


^o.»_2 g 


28.5±J;0 





Mn (GeV) 


0.94±0.28 


1 01+0-21 
-^•'^-^-0.21 


n 0-1-1-0.20 

U.»l_o.20 


mo 


Ma (GeV) 


1.12±0.34 


1 91-1-0.20 
-^•^-^-0.20 


90-1-0.20 

u.a»_o 20 


mo 


Ms (GeV) 


1.19±0.36 




1 90+0.20 
-^•^■^-0.20 


mo 


Ms (GeV) 


1.32±0.40 


1 34-1-0.21 

-^•■'*-0.21 


1 99+0-20 
-^•^^-0.20 


mo 



Table 5: Results by using the fits in Table |31 In the second column we give the values for the several observables 
as included in the fit. 
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Figure 1: From top to bottom and left to right, the ten panels correspond to the cross sections oi K p ^ 
K~p,K°n,'n-+T,~,Tv~T,+ ,Tr°T,°,TT°A, the 7r"S+ event distribution from K~p E+(1660)7r", the K~p -> 77A cross sec- 
tion, the 7r°E'' event distribution from the reaction K~p — >■ tt^tt^E" with px = 0.687 GeV and the total cross section of 
K^p — >■ tt^tt^E". The data points represented by black diamond, magenta square, orange circle, blue cross, cyan down- 
triangle and blue up-triangle in the first four panels are taken from Refs. [541 1471 [55l 1561 [57l 158] , respectively. The data in 
the fifth and sixth panels are from Ref. [59]. The 7r~E^ event distribution is taken from Ref. [24] and the K~p — >■ rjA cross 
section data are given in Ref. [25]. The measurements on the reaction K^p — > tt^tt^E^ are from Ref. [26]. The red solid lines 
and blue dashed lines represent the best fits from Fit I using Eqs. (|lip and p2l) (which is indicated by Fit I S), respectively. 
The areas covered by green hatched lines and the gray shaded areas correspond to our estimates of error bands for Fit I and 
Fit I S, in order. 
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Figure 2: The same as in Fig. [T] but for Fit II. We refer to Fig. [T] for notation. 
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Figure 3: The same as in Fig. [T]but for the C'(p)-Fit. We refer to Fig. [T]for notation. 
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^ 500 
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-250 
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Figure 4: Reversed sign energy siiift —AE and widtli of tiie Is state of Icaonic liydrogen. Tlie experimental data witli errors 
indicated by tiie sides of the rectangles correspond to the KEK [22], DEAR [21] and SIDDHARTA [23]. In the left panel, 
the empty square, circle and triangle are our results from Table[5]for Fit I, Fit II and C'(p)-Fit respectively, while the filled 
square, circle and triangle correspond to the results from Table[S]for Fit I, Fit II and C'(p)-Fit, in order. In the right panel, 
we show our final weight-averaged results in Table [7] by the (red) empty circle (denoted as OURS) and those from Ref. [16] 
(denoted as IHW in black plus symbol), Ref. [17] (denoted as CS in green triangle) and Ref. [T^ (denoted as MM in blue 



03 
O 



< 




5 6 7 

Different observables 



Figure 5: We show from left to right 7, Rc, Rn, flo+i ^tA, cr-n-iv and the baryon masses for A'^, A, E and H. For each 
observable the left most point corresponds to Fit I, next, going from left to right, one has the experimental value and finally 
the result for Fit II. For every of the fits we show the results obtained by employing each definition of x^, Eqs. (1111) (empty 
symbols) and H12p (filled symbols), as indicated in the figure. 
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1.38 1.39 1.4 1.41 1.42 1.43 1.44 1.45 1.38 1.39 1.4 1.41 1.42 1.43 1.44 1.45 

^ (McV) ^ (McV) 

Figure 6: Extrapolation of tlie K~p — ^ K~p amplitude T[y/s) in Eq. (|10|l to the subthreshold energy region for Fit I (red 
solid line) and Fit II (dashed blue line) in Table [T] The real part is shown in the left panel and the imaginary part in the 
right one. The shaded and hatched areas surrounding the solid and dashed lines, respectively, correspond to the statistical 
error bands from the fit in Table [1] The points with error bars correspond to the scattering length directly obtained from 
the SIDDHARTA measurement and the modified Deser-type formula, Eq. 



3.2 Results using the defined in Eq. ( [T2|) 

Next, we consider the stability of our fits if the is defined aceording to the standard definition Eq. instead 
of using Eq. ([TT]) . If a big change is observed by doing this modification in the outcome of a fit, then this clearly 
indicates that arbitrariness in the way chosen to describe nature is affecting our results. If this is the case such fit 
should be discredited. We conceive this check as a stability criterion. The fitted parameters for these new fits are 
shown in Table |4l At the semiquantitative level they are quite similar to those given in Table [TJ The only exception 
is fli for Fit I, although this parameter appears unstable when changing from one fit to the other within the same 
definition of x^- This is already the case in Table [TJ For the cross sections and event distributions shown in Figs.[Tl 

and 131 the new curves, given by the dashed lines, differ very little from the solid lines obtained previously. The 
uncertainty in the new curves is indicated by the shaded band around each dashed line. Then, there is stability 
under the change in the definition of the for these observables. We now consider those magnitudes given in 
Table [2] The new values employing the standard definition Eq. (I12|) are given in Table [5] and are shown pictorially 
in Figs. m and [5] by the filled symbols (squares, circles and triangles). Taking into account the error bars the results 
are compatible between the two definitions of the x^- The only exception happens for aj^_ in Fit II, for which the 
value obtained employing the definition of the x^ in Eq- (HH) is clearly different to the value obtained previously. 
Nevertheless, since the systematic uncertainty for our calculation of aj_^ is big, as indicated by the large error bar 
attached to the cross in Fig. |4j we consider that this deficiency is not significant. Regarding the C'(p)-Fit, not 
shown in Fig. [Sj one can observe by comparing the fifth columns in Tables [5] and [S] that the variation of all the 
quantities is within the one-sigma range. 

For the energy shift and width of the Is state of kaonic hydrogen we see clearly in the left panel of Fig U] that 
the new values, given by the filled square (Fit I), circle (Fit II) and triangle (C'(p)-Fit) are within the estimated 
errors compatible with those obtained before employing the definition Eq. (jlip . Nevertheless, the central values 
for each fit move around one sigma from each other when changing the definition of x^ • We also give in Table [6] 
the values of the scattering lengths a^-p, 07=0 and a/=i for these new fits employing the definition Eq. (jl2p for 
the x^- They are perfectly compatible to those given before in Table [3] In summary, our stability criterion is well 
fulfilled. 

From the discussion above, we take as resulting values from our study for the KN scattering lengths, /S.E and 
r the mean and variance calculated from the four NLO fits. That is. Fit I and Fit II in Tables [T] and E) For each 
quantity we add in quadrature the variance obtained and the largest of the statistical errors resulting from the x^ 
definitions Eq. (|12p and Eq. (Ilip . this gives our final estimate for the error bar. In this way, the spread in the 
values by changing the x^ definition is taken into account as a source of systematic uncertainty in our results given 
in Table [71 Our final values for A£' and F are also plotted in the right panel of Fig. [Hby the empty circle. There 
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Observable 


Fit I 


Fit II 


0(p)-Fit 




-0.75tro^+n.OOlH«3 


_Q 74+0.07 .-ooo+0. 09 
^•'^-0.08 ^ '^•='■'-0.08 


^^■"'-0.09 + ' -^■'^'-'-0.12 


07=0 (fm) 


-1.86t«;lt 1.351°;- 


_1 70+0.13 ng+0.18 


_1 72+0-18 +jl 47+0-30 
-^■'^-0.17 ^ '-^-^'-0.23 


07=1 (fm) 


0.37l°:°t + i0.65lH^ 


(, 01 +0.05 , ^Q+0.05 


0.381°:°? + i0.52in? 



Table 6: /S' p scattering length, aj^-p, and / = 0, 1 KN scattering lengths, a7=o, a7=i, in order, for Fit I, Fit II and 
C'(p)-Fit in Table H 



we also show some other recent determinations [151 EH EZl HI] that include SIDDHARTA data [23] in their fits. 



AE (eV) 


r (eV) 


"■K-p (fm) 


0.1=0 (fm) 


0.1=1 (fm) 


308 ± 56 


619 ± 73 


(-0.69 ± 0.16) + i (0.94 ± 0.11) 


(-1.74 ± 0.34) + i (1.31 ± 0.20) 


(0.36 ± 0.12) + i (0.56 ± 0.12) 



Table 7: Our final results for the quantities indicated in the first row. The error bars are calculated as explained in the 
text. 



4 Associated spectroscopy 

Now we analyze the baryon spectroscopy that stems from both fits I and II. We discuss only the results obtained 
from the NLO fits in Tabled] employing the definition Eq. (jlip . because those obtained from the NLO fits in Tabled 
are almost coincidentjf] To perform this study, we need to extrapolate the meson-baryon scatting amplitudes from 
the physical or first Riemann sheet to the unphysical Riemann sheets in the complex energy plane. The physical 
Riemann sheet is such that the imaginary part of the modulus of the three-momentum associated with every 
channel is positive. The other Riemann sheets are defined depending on which three-momenta are evaluated in 
the other sheet of the square root, with an additional minus sign. In practice the change of sheet can be easily 
performed by adding a new term to g{s)i, Eq. so that this function in its second Riemann sheet, g{s)iji, is 
given by [38] 



9{s)i + 
AttW 



ip{s) 



(14) 



with q{s)i-j the CM three-momentum of the ith state calculated in its first Riemann sheet, with positive imaginary 
part. Notice that in this way along the real axis and above the ith threshold the imaginary part of g{s)i-ji changes 
sign in the associated second Riemann sheet (denoted by II) with respect to first Riemann sheet. Let us denote 
the physical sheet as (-I-, . . .), where we have ten entries inside the bracket corresponding to the sign of 

the imaginary part of every qi{s) function for the ten coupled channels. Strictly speaking there are 10 unphysical 
Riemann sheets that are connected continuously with the physical sheet by crossing from s + iO"*" to s — iO^, which 
can be symbolized as (— , -f , . . .), (— , — , . . .), (—,—,—,-!-,...) and so on, adding an additional minus 
sign every time a threshold is passed over. Nevertheless, we notice that the gaps between the three thresholds 
TT^S", 7r~I]+ and tt+E" are quite narrow, and the same applies to the K^p, K^n and K^'EP , i^T+S" channels. 
Indeed, every of the sets of thresholds indicated would be degenerate if isospin symmetry were not broken in 
the pseudoscalar and baryon masses. In the following discussions, we do not distinguish the small differences of 
thresholds inside the three clusters. In such a way, we only have six unphysical Riemann sheets that are directly 
connected to the first one. We denote the second Riemann sheet as (— ,-|-,-|-,-|-,-|-,-|-,-|-,-|-,-|-,-l-), meaning that 
the sign of the imaginary part of qi is changed and becomes negative. The 3rd (3RS), 4th (4RS), 5th (5RS), 6th 
(6RS) and 7th (7RS) Riemann sheets correspond to (—,—,—,—, -f, -f, -f), (—,—,—,—,—,—,+,+,+,+), 



*We have explicitly checked that for the 0(p)-fit, except for the A(1670) resonance that is not generated, the rest of isoscalar and 
isovector resonances are quite similar to those from the NLO fits and no further insight results. 
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(-, -, -, -, -, -, -, +, +, +), (-,-,-,-,-,-,-,-,+, +) and (-,-,-,-,-,-,-,-,-,-), respectively. Apart 
from the pole position sji in the complex energy plane, the resonance R is also characterized by its couplings, which 
correspond to the residues (3i of the resonance pole and are calculated through 

T = _ lim . (15) 

s-^SR s — Sr 

We summarize the resonance pole positions and their couplings in Table [8] for Fit I and Table [9] for Fit II. 

The two relevant resonances with / = 0, namely A(1405) and A(1670), clearly show up in both fits The two poles 
corresponding to the A(1405) are found in the 3RS for both Fit I and Fit II. The narrower poles, i.e. those with 
small imaginary parts, from the two fits are perfectly consistent, while the broader ones differ slightly (taking into 
account errors) from the two fits. The narrower poles from our study also agree well with the recent determinations 
from Refs. [15l[T6l|8]. While for the broader poles, our results from Fit II are compatible with those in Refs. [151 [T6] . 
that also use the physical decay constants for tt, K and rj as in our Fit II. Concerning the broader pole, none of 
our results seem compatible with that in Ref. llBj , though we find that the real part of the broader pole in the case 
of Fit I is slightly above the K~p threshold. As pointed out in Ref. [12] the broader pole is much more dependent 
on the details of the theoretical approach than the narrower one. Our results for the broader pole are well inside 
the interval of values obtained in the same reference. Despite the broader pole changes its pole position according 
to the fit considered, the ttE and KN amplitudes have their peak values at nearly the same energy. This is due to 
the fact that the gradient of the modulus of the amplitudes is tilted so that it is not perpendicular to the real axis 
from the pole positions. This is clearly shown by the contour plots in Fig. [7]where the modulus of the elastic 7 = 
KN (blue dashed) and ttS (red solid) S-waves for Fit I (left panel) and Fit II (right panel) are plotted. For Fit 
I the gradient is tilted to the left from the broader pole position while for Fit II it is tilted to the right from the 
analogous pole, being always oriented towards the narrower pole. However, the gradient is always tilted to the left 
of the narrower pole. This figure also shows that for the ttT, state the two poles play a role in its shape on the real 
axis while KN is dominated by the narrower pole. In addition, the change in the position of the broader pole has 
more influence below the KN threshold and this is the origin of the differences between Fit I and Fit II in Fig. [HI 
as can be easily inferred from Fig. [7] considering the dashed lines. 




R,,l^ (GcV) 

Figure 7: Contour plot for the modulus of the elastic / = KN (blue dashed) and ttE (red solid) S-waves for Fit I (left 
panel) and Fit II (right panel). 

In addition to the pole positions, we also calculate the couplings of the resonances, see Tables [S] and O The 
first conclusion we can make for the A(1405) is that both the narrower and broader poles couple mostly to the ttS 
and KN channels. For the broader poles, the coupling strength to ttS is around twice as large as the one for the 
narrower pole, while the coupling to KN for both the narrower and broader poles is quite the same. Taking into 
account that the narrower pole gives strength around the KN threshold it dominates a reaction with production of 
KN, as pointed out above regarding Fig. [71 In turn, the signal for ttS is more dominated by the broader pole that 
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couples more strongly to this channel than the other narrower and higher in mass pole, see also Fig.[71 Concerning 
the channels with thresholds larger than the resonance masses, such as 77A and i^S, the narrower and broader poles 
couple to them weakly, much less than to ttE and KN . 

Next we consider the A(1670) resonance. The properties of this resonance are mainly determined by the K~p — >■ 
rjh. reaction that is included in our fits [3S], but not in the recent studies of Refs. [TSJ HH HZl [T^ . As as result these 
references do not obtain any information about the A(1670) resonance. In our case. Fit I and Fit II lead to quite 
similar good reproductions for the K~p — )■ TyA cross section, see the eighth panels in Figs. [1] and [21 respectively. 
As a result, the A(1670) pole positions and couplings from both fits perfectly agree between each other. They also 
agree with those of Ref. [8l[60] and with the PDG results [45]. Very similar poles in three different Riemann sheets, 
namely 3RS, 4RS and 5RS, are found. As explained in Ref. [HI [HI] all these poles reflect the same underlying 
resonance since they are connected continuously when changing the different Riemann sheets indicated in a soft 
way via a continuous parameter. We conclude that the A(1670) is most strongly coupled to the K'E. channel and 
has similar coupling strengths to tt'EjKN and 77 A channels. These couplings are also quite similar to those found 
in Ref. 8 . 

Regarding the I — 1 poles the resonance content is less clear since it depends on the details of the fits performed. 
For both NLO fits one has broad poles with central value masses of 1646 and 1686 MeV for fits I and II, respec- 
tively, that are in correspondence with the properties of the bumps associated with the E(1620) in the PDG [45] . 
Additionally, Fit II shows a pole at the central value position of 1741 — i 94 MeV with mass and width in the range 
given by the PDG [45] for the I](1750) resonance. From the residues shown in Table|9]its main decay channels are 
ttE, KN and rjT, as in [45]. Nonetheless, in our case the ryE decay width is smaller than the value suggested in [45] 
with large errors. The two last poles in Table IHI with 7=1 are connected in a continuous way by changing from 
the 5RS to the 6RS with a continuous parameter so that both poles correspond to the same physical object seen 
in different Riemann sheets [HUM]- In Fit II we also have narrow poles near the KN threshold. Similar poles were 
also reported in Refs. [51[S]. 

5 Conclusions 

We have studied the strangeness —1 5- wave meson-baryon scattering with 10 coupled channels by applying Uni- 
tary ChPT (a unitarization method based on an approximate algebraic solution to the N/D method) with the 
interaction kernel calculated up to 0{p^). The latter then corresponds to the NLO ChPT meson-baryon partial 
wave amplitudes. This study is prompted by the new precise measurement of the energy shift and width of the 
ground state of kaonic hydrogen. We have successfully reproduced a large amount of experimental data, which 
include the cross sections of K'p {K'p, K°n, tt+E", 7r-E+, 7r°E°, 7r°A, rjA, 7r°7r"i;°}, 7r~E+ and tt^E*' event 
distributions, three branching ratios in Eq. ([T]) measured at the K~p threshold, the ttA phase shift measured at the 
S~ mass and the new SIDDHARTA measurement on the energy shift and width of the kaonic hydrogen Is state 
Eq. We constrain further the fits by including, with rather large error bars, the pion-nucleon scalar isoscalar 
.S-wave scattering length, the nucleon a term and the masses of N, A, E and S. We confirm the consistency be- 
tween scattering data and the SIDDHARTA measurement, already highlighted in Refs. [TS] HH HZHIH] , although 
we consider additional recent data from the reactions K~p 77A [53] and K~p — >■ 7r°7r°E° [55], not included in 
those references. We have performed two types of fits. In Fit I we employ a common pseudoscalar decay constant, 
while in Fit II we distinguish the tt, K and ij weak decay constants. These two strategies are consistent with our 
0{p^) calculation of the interaction kernel since the difference between weak decay constants gives rise to higher 
orders in the counting. In addition we have studied for every fit the two definitions for the typically used in the 
literature for the study of strangeness —1 5- wave scattering. In one of them more weight is given to the presumably 
precise measurements of the energy shift and width of the Is state of kaonic hydrogen. We observe that the results 
are quite stable under the change of the definition for the x^- On the other hand, despite the knowledge of the 
K~p scattering length is an important ingredient to constrain the subthreshold extrapolation of the K~p S'-wave 
scattering amplitude, it is not enough to determine it in a precise manner. This is clear by the significant differ- 
ences in the subthreshold extrapolations resulting from Fit I and Fit II. In this sense, the systematic uncertainty 
affecting such extrapolation is larger than the statistical one reported in Refs. [TS] dH HZJ [TH] . The situation might 
be improved in the future once the interaction kernel is calculated at 0{p^) in Unitary ChPT. 

Discussions on the baryon resonance spectroscopy are carried out as well, for which we calculate both the 



17 



Pole 




/37rE|0 


l/3-.sli 


l/3.s|2 


l%jvlo 


l/3jfjvli 


I/3^a| 


l/3r,E| 


Ifeslo 


Ifesli 


A(1405) 
14361^(5 - * 1261^4 (3RS) 
14l7tl-i24tl (3RS) 


0+°-" 

'^■'-'-0.0 


s a+0.9 
»-»_0.4 


"■"-0.0 

1+°-° 

"■^-0.0 


0+°-° 

"■"-0.0 
"■"-0.0 


7 7+1.3 
'■'-0.7 
7 7+1.2 
'■ '-0.6 


"•"-0.0 
"•J^-0.0 


^•^-0.3 
-^•^-0.3 


"•"-0.0 
"•^-0.0 


9 1 +0.8 
^■-^-0.7 


"•"-0.0 

l+^'O 
"•^-0.0 


A(1670) 
16741^ -iStl 
1674j;| - i lltl {5RS) 
1673^1 -illtl (6RS) 


0+°-° 


o.8«:t 

q+°* 

U.9_g 2 

u.y_o 2 


0.0«:° 

"■"-0.0 


oo oo oo 
dd dd dd 
+ 1 + 1 + 1 
O O O 
odd 


1 t;+0.4 
^■^-0.2 

^■°-0.2 

^^■"-0.2 


"•"-0.0 

0+°-° 

"•"-0.0 

"•"-0.0 


-^•^-0.2 
1 7+0.5 
'-0.3 
-I 7+0.5 
-■-• '-0.3 


"•"-0.0 
"•"-0.0 
"•"-0.0 


10.8l°-J 

11 1+°-^ 

^^^■^-0.3 

-'--'-■-'--0.3 


l+"'° 
"•J^-0.0 

"•J^-0.0 


S / = 1 

imetliT - « 160t|^^ {4RS,5RS) 
18781^^ - i imtll (6RS) 


-,+1.4 
■^■^-0.5 

1 0+°'^ 


n n+°-° 
'J-^J-o.o 

0+°° 


■^■"-0.5 
c O+0.9 
^■»-0.6 


"■"-0.0 

0+°-° 

"■"-0.0 


"■"-0.0 
"■"-0.0 


^•^-0.3 
o 7+0.4 
■^•'-0.4 


"•"-0.0 
"•"-0.0 


'•^-1.2 
o g+1.1 
'-'•^-1.2 


"■"-0.0 

l+°^" 

"■^-0.0 





Table 8: Resonances from Fit I in Table[T] The positions of the resonance poles are in MeV and modulus of their couplings 
are given in GeV. The latter are represented by where the isospin label / is indicated when more than one isospin 

channel is possible. 
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Table 9: Resonances from Fit II in Table [T] For notation see Table [S] 



resonance pole positions and their residues (or resonance couplings). We find that the properties for the relevant 
isoscalar resonances are quite robust from the two fits. The two-pole structure for the A(1405) is confirmed, in 
which the properties of the narrower pole, close to the KN threshold, are quite stable. Our results are consistent 
with those of Refs. [TSl HH [HI For the other broader poles, somewhat small differences from the two fits are 
observed, with one pole above and the other below the K~p threshold. These broader poles depend more on the 
details of the theoretical approach and on the fit performed. Despite that, the peak positions on the real energy 
axis for the .^A'^ and ttS amplitudes are quite the same in both fits, because the different orientation of the gradient 
of the modulus of the amplitudes with respect to the broader pole position. The gradients always point towards 
the position of the narrower pole much closer to the energy axis. Nevertheless, a precise knowledge of the pole 
position of the broader pole is required in order to end with a precise extrapolation of the elastic K~p scattering 
amplitude below the KN threshold. The positions of the pole and couplings for the A(1670) from the two fits 
are nearly identical and in good agreement with the properties reported in the PDG. The isovector resonances 
depend more on the details of the fit. Nevertheless, in both cases one observes broad poles in the region of the 
E(1620) resonance. Fit II also reports a pole position that could be identified with the resonance E(1750). This 
fit accumulates strength around the KN threshold in / = 1 too, corresponding to two relatively narrow poles in 
this energy region. 
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